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Abstract 

Gauge effects on the fluctuation properties of the normal-to-superconducting 
phase transition in bulk and thin film superconductors are reviewed. Similar prob- 
lems in the description of other natural systems (liquid crystals, quantum field 
theory, early universe) are also discussed. The relatively strong gauge effects on the 
fluctuations of the ordering field at low spatial dimensionality D and, in particu- 
lar, in thin (quasi-2D) films are considered in details. A special attention is paid 
to the fluctuations of the gauge field. It is shown that the mechanism, in which 
these gauge fluctuations affect the phase transition order and other phase transition 
properties varies with the variation of spatial dimensionality D. The problem for 
the experimental confirmation of theoretical predictions about the order of phase 
transitions in gauge systems is discussed. Related topics: gauge effects on the criti- 
cal behavior of unconventional superconductors, disorder, quantum fluctuations in a 
close vicinity of ultra-low phase transition temperatures , are also briefly discussed. 



1 INTRODUCTION 

1.1 Gauge effects and superconductivity 

A remarkable example of gauge theory in condensed matter physics is the Ginzburg- 
Landau (GL) functional of superconductivity [HE]- The latter is invariant towards both 
global U(l) rotations of the order parameter field ip(x) and local gauge transformations 
(rotations) of the same field and the vector potential A(x) of magnetic induction B. 
These two properties of gauge invariance define the global and local U(l) symmetries of 
the GL theory. In both cases the gauge group is a one-dimensional Abelian continuous 
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group, i7(l). The spontaneous breaking of these global and local symmetries below the 
phase transition point allows the appearance of the superconducting phases: the uniform 
Meissner phase and, under certain circumstances, the mixed (Abrikosov vortex [3|) phase. 
While the Meissner phase [(0) 7^ 0, (B) = 0] is a mere product of the global {7(1) sym- 
metry breaking, the vortex phase [(i/j(x)) 7^ 0, (B(x)) 7^ 0] , where the equilibrium field 
configurations of both order parameter field ip(x) and magnetic induction B(x) are spa- 
tially nonuniform, occurs as a result of the spontaneous breaking of local {7(1) symmetry. 
The local gauge is important also for the description of magnetic field penetration in both 
Meissner and vortex phases. The penetration is described by an additional characteris- 
tic length, the London penetration length pfl [2], which is different from zero only when 
the symmetry is broken, i.e., in the ordered phases, where the equilibrium value of ip is 
different from zero. The local {7(1) symmetry is present when the superconductor is in 
external magnetic field, or, when this field is equal to zero but magnetic fluctuations exist 
and determine the vector potential A(x) as a purely fluctuating field. 

The fluctuations of relevant fields, ifi(x) and A(x) in usual superconductors are small and 
can be neglected. Therefore, the GL functional can be investigated in the lowest order 
mean field (MF) approximation (alias, "tree approximation" [HE]). This is the usual 
way of treatment of the GL free energy, in particular, the study of GL equations [TJ [2J, [3] . 
Within the tree approximation, the phase transition from normal to superconducting state 
in a zero magnetic field is of second order [U [2j 0] . 

For long time this phase transition has been considered as one of the best examples 
of second order phase transitions, which has an excellent description within MF. But 
in 1974 Halperin, Lubensky and Ma (HLM) [6] showed that the magnetic fluctuations 
change the order of the superconducting phase transition in a zero external magnetic 
field (Hq = \Hq\ = 0)), i.e., the order of the phase transition from normal-to- uniform 
(Meissner) superconducting state at T c0 = T C (H = 0) (see, also Ref. [7]). Since then 
this fluctuation change of the order of normal-to-superconducting phase transition (HLM 
effect) has been under debate. After 1974 an overwhelming amount of theoretical research 
on this topic has been performed but up to now there is no complete consensus about the 
order of this phase transition. 

In this review we shall consider some aspects of this problem but we shall not be able 
to discuss or mention all relevant contributions. However, it is important to emphasize 
that the investigation of the properties of the superconducting phase transition is also 
important to other areas of physics. Here we shall briefly enumerate and discuss several 
examples. 

1.2 Quantum field theory and other problems 

In elementary particle physics the gauge invariant theory similar to the GL theory of su- 
perconductivity is called Abelian Higgs Model jgl El UHl El US] ■ The same global {7(1) 
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and local U(l) gauge symmetries are present there but the phenomenon of spontaneous 
symmetry breaking occurs only for imaginary mass of the Higgs (ip— ) field, which is the 
analog of the field ip, describing the Cooper pair in a superconductor. In the absence of 
spontaneous symmetry breaking this model would describe an ordinary electrodynamics 
of charged scalars, but the situation becomes more interesting when the mentioned mass 
is imaginary and the breaking of symmetry occurs. Now the symmetry breaking phenom- 
ena receive other names and another physical interpretation. The breaking of both global 
and local gauge symmetries ensures a mechanism for a transformation of the two initial 
scalar fields - analogs of the components ip' and ip" of the complex field ip = ip' + and 
two massless photon fields - analogs of the two independent components Aj of the vector 
potential A = {A,;j = 1,2,3; V. A = 0} in a superconductor, to four massive particle 
fields: the so-called Higgs boson, which is analog of the spontaneous order > in a 
superconductor, and three massive vector field components, i.e. a massive three dimen- 
sional vector field. The mass of this new vector field is proportional to the electric charge 
\q\ and magnitude \ip\ of the Higgs field and, as a matter of fact, this is exactly the way, in 
which the London penetration length in a superconductor depends on the electron charge 
|e| and the modulus j^l of the superconducting order parameter ip. Thus the spontaneous 
breaking of the local gauge symmetry leads to the formation of massive particles without 
spoiling the gauge invariance of the theory. This is called "the Higgs mechanism" . It plays 
a fundamental role in the unified theory of electromagnetic, weak and strong interactions 
(see, e.g., Ref. [El [13]). 

It is easy to see that there is something quite common between the phenomena of sponta- 
neous breaking of the continuous symmetries in superconductivity theory and in quantum 
field theory. The superconducting phase > is the exact analog of the Higgs boson in 
the Abelian-Higgs model, whereas the appearance of a massive vector field has its analog 
in the finite London penetration length, mentioned above. Of course, there is no obsta- 
cles in interpreting the phenomena of spontaneous symmetry breaking in quantum field 
theory as phase transitions by taking the Higgs mode mass as tuning parameter. The 
phase transition will occur at zero mass of the Higgs boson. 

A similar phenomenon of spontaneous breaking of both global U(l) and local gauge sym- 
metries is possible also within the scalar electrodynamics due to mass insertions from the 
radiation corrections, as shown in Ref. [Hj. The radiation corrections, in analogy with 
the magnetic fluctuations in a superconductor, generate an imaginary mass to the initially 
massless scalar field in this theory and the latter becomes very similar to the Abelian- 
Higgs model. Here the symmetry breaking leads to the appearance of massive scalar and 
vector fields describing neutral scalar meson and vector meson, respectively [14J. The ra- 
diation corrections to the Lagrangian of the massless scalar electrodynamics [Hj resemble 
very much, in particular in their mathematical form at D — 4, the magnetic fluctuation 
corrections to the GL free energy of 4D superconductors [7j. The interrelationships be- 
tween the superconductivity theory and the gauge theories of elementary particles has 
been comprehensively discussed in Ref. |15j . Note also the relation between the supercon- 



3 



ductivity GL functional and the CP N 1 confinement model (see Ref. [TB]) and extensions 
to non-Abelian theories [PT] . 

The gauge theories, mentioned so far, and their extensions have a wide application in the 
description of the Early Universe [T5l [18]. Another interesting gauge theory is that of 
the nematic-to-smectic A phase transition in liquid crystals. According to the Kobayashi- 
McMillan-de Gennes theory [191 EH El] the smectic-74 order is described by two order 
parameters: the nematic director vector and the complex scalar which gives the long 
molecule centers of masses. When a description, analogous to that of superconductors, is 
introduced, as suggested by de Gennes [21] , the director vector is substituted with a gauge 
vector field, which is quite similar to the vector potential A in the GL functional. Apart 
from some specific features intended to take into account the liquid crystal anisotropy, 
the effective free energy of the smectic A looks very much like the GL free energy of 
superconductors (see, also, Refs. [22j (23J, El] ) . A gauge theory in condensed matter physics 
that has some (but not very close) similarity with the GL free energy of superconductors 
is the Chern-Simons-Ginzburg-Landau (CSGL) effective functional of the quantum Hall 
liquid state, in particular, in the description of phase transitions between the plateaus in 
the quantum Hall effect [25], EH1 EH EH] • Let us mention also the liquid metallic hydrogen, 
where the problem of the superconducting-to-superfluid phase transition [29] is also related 
to the topics, discussed in this review. 

1.3 About the investigation of the fluctuation-driven first order 
phase transition in superconductors 

Now we will focus on how the fluctuations effect the properties of the normal-to-superconducting 
phase transition in a zero external magnetic field, which is the subject of the present re- 
view and in particular cases we shall refer to related topics in other natural systems. We 
shall consider in more details the HLM effect in 2D and quasi-2D superconductors with 
a special emphasis on the problem of the theoretical predictions reliability. Following 
Ref. [6] we will use two theoretical methods. 

In Sec. II we use the MF like approximation of Ref. [6], where the spatial fluctuations 
of superconducting order parameter if) are neglected, and the magnetic fluctuations are 
taken into consideration. This treatment was justified for well established type I super- 
conductors, where the London penetration length A is much smaller than the coherence 
(correlation) length £. Due to the neglecting of spatial fluctuations of if>, the results of such 
MF treatment should be valid outside the Ginzburg critical region [4] . A weakly first order 
phase transition occurs as a result of new small |^| 3 -term, which appears in the effective 
free energy of 3D superconductors within the framework of the MF like approximation [5] . 
But this HLM effect was found to be very small and experimentally unobservable even 
for well established type I bulk (3D-) superconductors, such as Al, where the GL number 
k = A/£ = 10 -2 < 1. It has been recently shown §U\ M, M, G31 EU ESI [36] that the 
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HLM effect is much stronger in quasi-2D superconductors than in bulk (3D) samples. 
Moreover, as shown in this series of papers, the effect appears by a term of type | , 0| 2 hi|?/'| 
in the effective free energy, and essentially depends on the thickness L of quasi-2D su- 
perconducting films. These circumstances provide a real opportunity for an experimental 
verification of the effect - a topic of discussion throughout the present review, and in 
particular, in Sec. II. 

We must emphasize that MF results are not valid for thermodynamic states in the 
Ginzburg critical region (5T)q = \T G — T c0 \/T c0 around the equilibrium phase transi- 
tion point T c q. In certain classes of high-temperature superconductors the Ginzburg re- 
gion exceeds 0.1 K whereas in usual low-temperature superconductors it is very narrow, 
(ST)g ~ 1CT 12 — 10~ 16 K, and the respective critical states are experimentally unaccessi- 
ble (see, e.g., jlj). If the metastability states, which go along with the weakly first order 
phase transition predicted by MF, extend over temperature intervals larger than the size 
of the Ginzburg region, one can conclude that the MF prediction of the HLM effect is 
reliable. In Sec. II we justify the MF analysis reliability for element superconductors as 
Al, W, In. We will show that the condition for enough wideness of metastability regions 
is quite strong and, perhaps, irrelevant in experiments. We will also demonstrate that the 
magnitude of the critical magnetic field is crucial for HLM effect to be observed in thin 
superconducting films (see also the discussion in Sec. II. 9). 

In Sec. Ill we review some results for the HLM effect obtained with the help of the 
renormalization group (RG) method [HE]- The latter allows simultaneous treatment of 
both superconducting and magnetic fluctuations in the asymptotically close vicinity of 
phase transition point. The lack of fixed point of one- loop RG equations for conven- 
tional superconductors in a zero external magnetic field was interpreted as a signal for a 
fluctuation-driven first order phase transition [6]. This result shows that the local gauge 
magnetic fluctuations are relevant also in the Ginzburg region of strong -^-fluctuations, 
and under certain circumstances, they can change the phase transition order to a weakly 
first order, as is outside the critical region. The RG investigations of the superconducting 
phase transition order in a zero magnetic field has been recently reviewed in Ref. [37] . 
While the latter review emphasizes the investigation of the HLM effect in high orders 
of the loop expansion [U [5], here we shall stress on RG results concerning the effects 
of anisotropy, quenched disorder and extended symmetries of the Higgs ^-field on the 
phase transition properties. The RG study of unconventional superconductors will also 
be discussed (Sec. III. 7). In the remainder of this paper we will consider mainly the phase 
transition to the superconducting state but results for other systems and related topics 
will be also mentioned in brief; see, e.g., Sec. III. 3 - III. 6. 



5 



2 MF STUDIES 



2.1 General GL functional 



The GL free energy 
given in the form 



of D-dimensional superconductor of volume Vd = (Li...Ld) is 
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In Eq. (1) the first Landau parameter a = a^{T — T c q) is expressed by the critical tem- 
perature T c0 = T C (H = 0) in a zero external magnetic field (H = \H\), b > is the second 
Landau parameter and e = |e| is the electron charge. The square B 2 of the magnetic 
induction B = (H + AttM), is given by the vector potential A(x) = {Aj(x), j = 1, D} 
in the following way 

i A f().\j <).\, 1 

dxi 



B' 



D 

h J = 1 



(2) 



here the vector potential obeys the Coulomb gauge V ■ A(x) = 0. For 3D supercon- 
ductor the relation B = V x A(x) can be used and when B = B is uniform along the z- 
axis, the Landau gauge A (x) = B (—y/2,—x/2,0) can be applied. This representation 
can be generalized for D > 2 - dimensional systems, where the magnetic induction Bq is 
a second rank tensor l38l: 



Boij = B (5n5j 2 — 5j 2 5, 



a j 



(3) 



If we use the notation x = (xi,X2,r), where r is a (D — 2) - dimensional vector, perpen- 
dicular to the plane (xi, x 2 ), in the 3D case we have f— (0, 0, z), and 



-^i = -jtjkiBoki = B Q 5 j3 , (4) 

where Ejki is the antisymmetric Levi-Civita symbol. The Landau gauge and Eqs. (J3j) - (SJ) 
can be used for uniform B = Bo when 6B - fluctuations are neglected. In the prevailing 
part of our study we will apply the general Coulomb gauge of the field A(x), which does 
not exclude spatial dependent magnetic fluctuations 5B(x). 

In nonmagnetic superconductors, where the mean value (M) = (M — 5M) of magnetiza- 
tion M is equal to zero in the normal state in a zero external magnetic field, the magnetic 
induction in presence of external magnetic field takes the form: 

B = H + 8H(x) + 4tt5M(x) , (5) 

where H is the (uniform) regular part of the external magnetic field and 5H is an irregular 
part of H created by uncontrollable effects. We neglect the irregular part 5H and set 
H = 0, then B contains only a fluctuation part B = 5B{x) = 4irdM(x) that describes 



6 



the diamagnetic variations of M(x) around the zero value (M) = due to fluctuations 
6ip(x) of the ordering field ip(x) above (T > T c0 ) and below (T < T c0 ) the normal-to- 
superconducting transition at T c0 . Note, that the non-fluctuation part A = [A(x)—SA(x)] 
corresponds to the regular part B = (H +(M)) = of B in nonmagnetic superconductors 
((M) = 0) in a zero external magnetic field (Hq = 0). Then we can set Aq(x) = and, 
hence, 5A(x) = A(x), so we have an entirely fluctuation vector potential A(x), which 
interacts with the order parameter ip(x). This interaction can be of type |?/'| 2 ^4 and 
|-^| 2 y4 2 and generates all effects discussed in the paper. 

We accept periodic boundary conditions for the superconductor surface. This means 
to ignore the surface energy including the additional energy due to the magnetic field 
penetration in the surface layer of thickness equal to the London penetration depth A(T) = 
Ao|^o|~ 1//2 > = \T — T c0 |/T c0 ; A = (mc 2 b/87ie 2 a T c0 ) 1 ^' 2 is the "zero-temperature" value 
of A. This approximation is adequate for superconductors of thickness L ^> A(T) 
a , where a is the lattice constant and L = min{Li, i = 1, ...,£)}. As we suppose 
the external magnetic field to be zero (Hq = 0) or very small in real experiments, the 
requirement L ^> A(T) can be ignored and we have the simple condition L ^> a . 

In microscopic models of periodic structures the periodic boundary conditions confine the 
wave vectors k^ = {ki = (2irni/Li); i = 1,...,D} in the first Brillouin zone [— (n/cio) < 
ki < (7r/ao)] and the expansion of their values beyond this zone can be made either by 
neglecting the periodicity of the crystal structure or on the basis of assumption that large 
wave numbers k = \k\ have a negligible contribution to the calculated quantities. The last 
argument is widely accepted in the phase transition theory, where the long-wavelength 
limit (ka 1) can be used. In particular, this argument is valid in the continuum 
limit (Vd/clo — > oo). Therefore, for both crystal and nonperiodic structures we can use 
the cutoff A ~ (vr/ao) and afterwards extend this cutoff to infinity, provided the main 
contributions in the summations over k come from the relatively small wave numbers 
(k A). In fact, this is a quasi-macroscopic description based on the GL functional (1), 
so the microscopic phenomena are excluded from our consideration. 

The GL free energy functional takes into account phenomena with characteristic lengths 
£o and Ao or larger (£ and A), where A(T) is the London penetration length, mentioned 
above, and £(T) = £oN -1 ^ 2 is the coherence length [2]; here £ — (^ 2 /4ma T c0 ) 1//2 is the 
zero-temperature coherence length. In low-temperature superconductors £o and Ao are 
much bigger than the lattice constant ao- Having in mind this argument we will assume 
in our investigation that A (jr/ao). Whether the upper cutoff A is chosen to be either 
A ~ 1/£q or A ~ 1/Ao is a problem that has to be solved by additional considerations. 
According to arguments presented in Ref. [30] and Sec. II. 6, we will often make the choice 
A ~ Zq 1 . 

We will use the Fourier expansion 

M*) = -LYsMk)^* (6) 

Vd k 
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and 

"D k 

where the Fourier amplitudes Aj(k) obey the relation A* (k) = Aj(—k) and k.A(k) = 0. 
The Fourier amplitude ip(k) is not equal to ip*(—k) because if)(x) is a complex function. 
For the same reason ip(Q) = ip(k = 0) is a complex number. 

The functional (1) is invariant under global £7(1) rotations, defined by ij){x) — > T/>(x)exp(ia), 
where the angle a does not depend on the spatial vector x, and under the local U(l) gauge 
transformations if)(x) — > ip(x)exp[ia(x)], A(x) — > A(x) + (hc/2e)'Va(x). According to the 
discussion in Sec. I, we have to investigate the spontaneous breaking of these symmetries, 
that is, the ordered phases and the phase transitions in the superconductor. The HLM 
effect, on which we are going to focus the attention, is one of the results of the local U(l) 
symmetry breaking. 



2.2 Notes about the MF like approximation 

The effect of the superconducting fluctuations 8ip{x) on the phase transition properties 
is very weak, as in usual superconductors, and is restricted in a negligibly small vicinity 
(|t | ~ 1CT 12 -T- 1CT 16 ) of temperature T c0 , we will assume that Sip(x) = 0, i.e., ip ~ 
(ijj(x))] from now on we will denote (ip(x)) by ip. Therefore, we apply the mean-field 
approximation with respect to the order parameter ip(x). Within this approximation we 
will take into account the £y4(x)-fluctuations of Bq = 0, i.e., A{x) = 5A(x). Furthermore, 
the y4(x)-fluctuations can be integrated out from the partition function, defined by: 

Z{i>) = J VAe' F ^ A)/kBT , (8) 

where the functional integral J T>A is given by 

/oo D 
n n dAjwsfrvA*)) ■ (9) 
00 j=l x&V D 

The integration is over all possible configurations of the field A(x); the 5-function takes 
into account the Coulomb gauge. 

The partition function Z(ip) corresponds to an effective free energy Tt> of the .D- dimensional 
system: 

F D = -k B T\nZ(ip). (10) 

The magnetic fluctuations will be completely taken into account, if only we are able 
to solve exactly the integral (181) . The exact solution can be done for a uniform order 
parameter ip. The uniform value of ip is different from the mean- field value of i/j, because 
the uniform fluctuations of ip{x) always exist, so we should choose one of these two 
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possibilities [30,, 35]. The problem of this choice arises after calculating the integral (jHD 
at the next stage of consideration, when the effective free energy Tx> is analyzed and 
the properties of the superconducting phase (if) > 0) are investigated. The effective free 
energy is a particular case of the effective thermodynamic potential in the phase transition 
theory [H [5] and we must treat the uniform if> in the way prescribed in the field theory 
of phase transitions. It will become obvious from the next discussion that we will use a 
loop-like expansion, which can be exactly summed up to give a logarithmic dependence 
on \ip\. 

Due to the spontaneous symmetry breaking of the global U(l) continuous symmetry of 
the ground state if) ^ 0, the effective free energies discussed in this Section depend on 
the modulus of the complex number if) = \ip\e ie but not on the phase angle 9, which 
remains arbitrary. That is why we will consider the modulus as an "effective order 
parameter" as the angle 9 does not play any role in the phenomena investigated in this 
Section. The quantity \ip\ remains undetermined up to the stage when we determine 
the equilibrium order parameter \ip \ by the equation of state [OJ-'d (VO /9ip] = 0. This 
equation gives the equilibrium value ip of if) and the difference Sif) — ("00 - VO can be 
treated as the uniform (zero dimensional) fluctuation of field ip(x). The x-dependent 
fluctuations 8if>(x) have been neglected because of the if> uniformity. The solution ip will 
be stable towards the uniform fluctuation 5ip, provided the same solution if) = \ip \e ld ° 
corresponds to a stable (normal or superconducting) phase; the phase angle 9q remains 
unspecified. 

We begin our investigation by setting ip uniform but at some stage we will also ignore the 
uniform fluctuation Sif) and deal only with the equilibrium value if>o of if). The equilibrium 
value will be calculated after taking into account magnetic fluctuations, so it will be differ- 
ent from the usual result |^o| = (|a|/6) 1//2 |2J, where both magnetic and superconducting 
fluctuations are ignored. This simplest approximation for the equilibrium value of if) is 
obtained from the GL free energy (1), provided e = and the gradient term is neglected. 
Hereafter we will keep the symbol \ipo\ for the equilibrium order parameter in the more 
general case, when the magnetic fluctuations are not neglected, and will denote the same 
quantity for e = by 77 = \ipo(e = 0)| = (\a\/b) l l 2 . 

The above described approximation neglects the saddle point solutions of GL equations, 
where (if>(x)) is x-dependent. Therefore, the vortex state that is stable in type II super- 
conductors cannot be achieved. This is consistent with setting the external magnetic field 
to zero, so the vortex state cannot occur in any type superconductor. Our arguments can 
be easily verified with the help of GL equations [2] for a zero external magnetic field; the 
only nonzero solution for ip in this case is given by 77 = ( | a. | /6) 1//2 , although the magnetic 
fluctuations A(x) = 8A(x) are properly considered. 

In conclusion we can argue that the described method will be convenient for both type I 
and type II superconductors in a zero external magnetic field, if the ^-fluctuations have a 
negligibly small effect on phase transition properties T c0 = T C (H Q = 0), where T c denotes 
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the phase transition line for any H > 0. For type II superconductors in H > 0, two 
lines T c i(H ) and T c2 (H ) should be defined, usually given by H cl (T) and H c2 {T) [2]. 

2.3 Effective free energy 

When the order parameter ip is uniform the functional (1) is reduced to 

F(ij>,A)=F (il>) + F A (iJ>) (11) 

with 



F ^) = V D (a\^\ 2 + ^) (12) 



and 



';3 



F aW = -jd- X { p{^)A\x) + 1 ,yJ^- d 0> ■ (13) 



Here p = po\ip\ and po = (Bwe /mc ). It is convenient to calculate the partition function 
Z(ip) and the effective free energy in the fc-space, where Eqs. ([9]) and ( |T3l) take the 
form 

D k<A 

/ FI II dReAjfydlmAjfyS k ■ A(k) (14) 

F A {^) = F A (0) + AF A &) . (15) 

F A (0) = ^J2 k2 \Mk) 2 , (16) 

af a &) = P Y,\Mk)\ 2 ; (17) 



and 



Here 



and 



note, that we have used the Coulomb gauge k.A(k) = 0. 
Then the partition function (jSJ) becomes 

Z(iP) = e- FoW/kBT Z A (^j) , (18) 

where 

Z A (ip) = J VAe~ FAW/kBT (19) 

with F A (ip) given by (TT5"|) and the functional integration defined by the rule (I14p . With 
the help of Eqs. ffTUl) - f|T9|) the effective free energy Foi.^) will be 

^(V0 = ^o(V0+^>(V0> (20) 



10 



where F (ip) is given by Eq. (jT2"|) and 

Tfty) = -k B T\n 



Z(0) 



(21) 



is the ■0- dependent fluctuation part of J- (if)). In Eq. (1201) the -^-independent fluctuation 
energy {—k B T\n [^(0)]} has been omitted. This energy should be ascribed to the normal 
state of the superconductor, which by convention is set equal to zero. 



Defining the statistical averages as 

((•••)) 

we can write Eq. (I2ip in the form 



jVAe~ FA{ W kBT {...) 



T f (if)) = -k B T\n ( e -^A(i>)/k B Ty 



(22) 



(23) 



Eq. (|23|) is a good starting point for the perturbation calculation of !Ff(ifj). We expand 
the exponent in Eq. (|23|) and also take into account the effect of the logarithm on the 
infinite series. In result we obtain 



00 C-1V 



(24) 



where (...) c denotes connected averages [1]. Now we have to calculate averages of type 



(A Q (/c 1 ),A /3 (fc 2 )...A 7 (fc n )) c . 
Here we will use the Wick theorem and the correlation function of form 



c\f\k,k') = (Mk)M-k')) = SrpGfm 



where 



GfJk) = (Ai(k)Aj(-k)) 



Ank B T 



(25) 

(26) 
(27) 



and ki = (ki/k). 

The calculation of the lowest order terms (I = 1,2,3) in Eq. (|24p with the help of 
( 12"T|) is straightforward. The perturbation terms in (12~4"|) are shown by diagrams in Fig. 1. 
The infinite series ( |24|) can be exactly summed up and the result is the following logarith- 
mic function 



i + 



k 2 



[28) 



The same result for Fffy) can be obtained by direct calculation of the Gaussian functional 
integral ([8]). This is done using the integral representation of 5-function in (Q or (Tfj 
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Figure 1: Diagrammatic representation of the series (24); • represents the p- vertex in (13) 
and (17), and the solid lines represent bare correlation functions (|A,-(A;)| 2 ). 

but it introduces an additional functional integration that should be carried out after the 
integration over Aj(x). 

Eqs. ( flOl) . ( |20|) and ( 1281) give the effective free energy density 

f D ^)=FD^)/V D (29) 

in the form 

ZdW = /oW + A/ d (V0, (30) 

where 

/o(0 = «W + ^| 4 (31) 

and 

A/D( , H (£_^ Eln(l + iL). (32) 

k 

Eqs. (1201) and f[2"9~j) - (13"2"j) are the basis of our further considerations. We should mention 
that the fluctuation contribution A/d^) to f{ip) transforms to convergent integral in the 
continuum limit 




where Kb = 2 1 D ir D ^ 2 /T(D/2) for all spatial dimensionalities D > 2. But the terms 
in the expansion of the logarithm in (13"2l are power-type divergent with the exception 
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of several low-order terms in certain dimensionalities D. Therefore, we will work with 
a finite sum of an infinite series with infinite terms. In our further calculations we will 
keep the cutoff A finite for all relevant terms in Afjj(ij}). This is the condition to obtain 
correct results. 

2.4 2D-3D crossover 

The dimensional (2D-3D) crossover has been considered in Refs. [34"l 139] . Here we follow 
Ref. [34J, where the effective free energy density f(if>) = fzi^j) = -^(VO/^ of a thin 
superconducting film of thickness L and volume V = V 3 = (L Q LiL 2 ) is derived in a more 
general way, which allows the investigation of the dimensional crossover. Now one should 
perform the integration (33) only with respect to the wave vector components k\ and 
&2, corresponding to the large sizes Lj ^> L (J = 1,2). The result for the effective free 
energy density is [34] : 

m 

where 

is given by the sum 

S 

and q = \q\, q= (h, k 2 ). 

In Eqs. (J34l) -(j36l). the integral J(p) and the sum S(q,p) over the wave vector k = (q, k ) 
are truncated by the upper cutoffs A and A . The finite cutoff A is introduced for the 
wave number q and A stands for fc . 

As our study is based on the quasimacroscopic GL approach the second cutoff A should 
be again related to £o rather than to the lattice constant ao, i.e. Ao ~ (l/£o), which means 
that phenomena at distances shorter than £o are excluded from our consideration. We 
will assume that the lowest possible value of A is (vr/^ ), as is for A, but we will keep in 
mind that both A and A can be extended to infinity, provided the main contribution to 
the integral J(p) and the sum 5* comes from the long wavelength limit (g£o <C 1). 

In a close vicinity of the phase transition point T c q from normal {ip = 0) to Meissner state 
(|^| > 0) the parameter p ~ \ip\ 2 is small and the main contribution to the free energy f(ijj) 
will be given by the terms in S with small wave vectors k -C A. This allows an approximate 
but reliable treatment of the 2D-3D crossover by expanding the summation over /c in (1361) 
to infinity - Ao ~ oo. A variant of the theory when Ao is kept finite (A = Ao = 7r/£o) can 
also be developed but the results are too complicated [36]. Performing the summation 



a\^\ 2 + -\^ + k B TJ[p(^)} 



(34) 



J(p)= I t^qS (q,p) 



(35) 



+Ao 

- Y In 

fc =-A 



1 + 



q 2 + k 2 



(36) 
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and the integration in Eqs. (I55]) - (I35|) we obtain J(p) = (A 2 /27rL )/(p), where 



Hp) 



dy In 



sh {\ L Q k^p~Ty) 
sh (| L A^y) 



(37) 



The integral ( 1371) has a logarithmic divergence that corresponds to the infinite contribution 
of magnetic fluctuations to the free energy of normal phase (T c0 > 0, </? = 0). Such type 
of divergence is a common property of many phase transition models. In the present 
case, as is in other systems, this divergence is irrelevant, because the divergent term 
does not depend on the order parameter ip and the free energy f(ip) is defined as the 
difference between the total free energies of the superconducting and normal phases: 
fW = (fs-f N ). 

Introducing a dimensionless order parameter ip = (ip/?p ), where ipo — (ci T c0 /bY^ 2 is the 
value of ip at T — 0, we obtain the free energy (134"|) in the form 



/(¥>) 



H c0 

8tt 



2V + -^| 4 + 2(l + t )C% 5 



(3* 



with I([j,<p 2 ), given by Eq. (1371) . \i = (1/hk,) 2 , A = tt/£ , and 

2n 2 k B T c0 



C 



Lo^H 2 



(39) 



From the equation of state (df /dip = 0) we find two possible phases: (p 00 = and the 
superconducting phase (ip > 0), defined by the equation 



. 2,(1 + ^o)C£o£o , 2 \ 
4ttA 







where 



K{z) 



1 dy — h ^ L o A VlF^) 



(40) 



(41) 



o 



The analysis of the stability condition (d 2 f/dip 2 > 0) shows that the normal phase is a 
minimum of f(<p) for to > 0, whereas the superconducting phase is a minimum of /(<£>), if 



1 > ^(l+t )CL Afi 2 K(ij ri 



2 ) 



where 



K(z) 



dy 
o y + z 



coth (| L Ay/y + 



+ 



2sh 2 (| L A^y-+^) 



(42) 



(43) 



The entropy jump AS" = (AS/V) = [—df((po)/dT] per unit volume at the equilibrium 
point of phase transition T c ^ T c0 is obtained in the form 



AsS{T c 



1 + 



Clip 



cQ) 



(44) 
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where (p c $ = <po(T c ) is the jump of the dimensionless order parameter at T c . 

The second term in AS can be neglected. In fact, taking into account the equation 
f[(p (T c )] = for the equilibrium phase transition point T c we obtain that |C/(<p )/Vol is 
approximately equal to \t cQ +tp 2 . /2\, where ip 2 cQ and the dimensionless shift of the transition 
temperature t c o = to(T c ) are expected to be much smaller than unity. The latent heat 
Q = T C AS(T C ) and the jump of the specific heat capacity at T c , AC = T c (dAS/dT) can 
be easily calculated with the help of Eq. (j44p . For this purpose we need the function 
(fo(T), which cannot be obtained in analytical form from Eq. fj4"0]) . 

Eqs. (38) and (40) can be analyzed numerically. These relatively simple 2D-3D crossover 
formulae can be used in investigations of specific substances by variation of thickness Lq of 
films from Lq ^> A -1 (3D system) to ao < Lq <C A -1 (quasi-2D system) and even to a 2D 
system for L = a Q . Then one can vary the effective dimension of the system D e ff(L A) 
as a function of L Q A Q from D = 2 to D = 3 [43J. However, from a purely calculational 
point of view we have found more convenient to consider particular dimensions of interest 
separately and then to compare the results in order to demonstrate the relevant differences 
between the bulk (3-D) and thin film properties. This approach is applied below. 



2.5 Effective free energy for particular dimensions 

For purely 2D superconductor, consisting of a single atomic layer, we can use Eqs. ff29|) - 
( l32l) by setting D = 2, and calculate Af 2 (ip) with the help of rule ( |33l) : 



k B T 



57T 



Po\^?\ ,,,121 ( Po\^? 



(A^ + p h/f)ln 1 + ^- -poH 2 hi 



A 2 J ruiri V A 2 
The first term of the above free energy can be expanded in powers of \ip\ 2 : 



(45) 



AA W = (M) {^ + ^.n(-^) + ^ . (46) 



Thus we obtain the result from Ref. [40]. This case is of special interest because of the 
logarithmic term in the Landau expansion for f(ip), but it has no practical application 
for the lack of ordering in purely 2D superconductors. 

For quasi-2D superconductors we assume that (27r/A) > L ^> a , where L is the 
thickness of the superconducting film, and the more precise choice of the upper cutoff 
A ^ (1/ a ) for the wave numbers fcj is a matter of additional investigation [30]. In order 
to justify this definition of quasi-2D system one can use the 2D-3D crossover description 
presented in Sec. 2.4. The summation over the wave number k = (2nn /L ) in Eq. (36) 
cannot be substituted with an integration because L <C Lj and the dimension L does 
not obey the conditions, valid for Lj, (J = 1,2) [H], 021 H3]. Therefore, for such 3D 
system we must sum over fc an d integrate over two other components (k\ and k?) of the 
wave vector k (see Sec. 2.4). This gives an opportunity for a systematic description of 
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the 2D-3D crossover as shown in Sec. 2.4. In the limiting case of very small thickness 
the 2D-3D crossover theory (Sec. 2.4) leads to a result, which is obtained more simply 
in an alternative way, namely, by ignoring all terms corresponding to ho 7^ in the sum 
in Eq. ( |32j) . This corresponds to the supposition that the quasi-2D film thickness can- 
not exceed 2nA. Assuming this point of view the real physical size of the quasi-2D film 
thickness will depend on the choice of the cutoff A. It is certain at this stage that A > £ 0) 
because the (quasi)phenomenological GL theory does not account phenomena of size less 
than £o- The upper cutoff A of wave numbers can be defined in a more precise way at 
next stages of consideration; see Sec. II. 6 - II. 8. 

For a quasi-2D film we have the expression: 

A/(V0 = ^-A/ 2 (V0 , (47) 

where A/ 2 (V0 is given by Eq. fT45|) . 

For bulk (3D) superconductor we obtain: 



A/ a (V0 



2tt 



Po|# 



• (48) 



The Landau expansion in powers of in this form of fs(ip) confirms the respective 
results of Refs. [6j[7], moreover it correctly gives the term of type Pq|?/>| 4 , which has been 
considered small and neglected in these papers. 

For 4D-systems A/d(V0 becomes 



A>,| ( f' + A'ln(l + 4*V^| 4 lnfl+ ' 



A 2 J rUir ' V PoM 



(49) 



The above expression for A/4 (ip) can be also expanded in powers of to show that 
it contains a term of the type \ip\ 4 In (a/PoI^I/A), which produces a first order phase 
transition; this case is considered in the scalar electrodynamics [H], as mentioned in 
Sec. I. In our further investigation we will focus our attention on 3D and quasi-2D 
superconductors. 

The free energy density A/d(V0 can be expanded in powers of \ip\ but the Landau ex- 
pansion can be done only in an incomplete way for even spatial dimensions. Thus ^(VOj 
f^ip), and f(ip) - the free energy density corresponding to the quasi-2D films, contain 
logarithmic terms, which should be kept in their original form in the further treatment of 
the function A/d(^/>) in the Landau expansion. The analysis has been performed [35] in 
two ways: with and without Landau expansion of A/d^). These variants of the theory 
are called "exact" theory (ET) and "Landau" theory (LT), respectively [33J. It has been 
shown [35] that these two ways of investigation give the same results in all cases except 
for quasi-2D films with relatively small thicknesses (Lq <C £o)- It seems important to es- 
tablish the differences between two variants of the theory because the HLM effect is very 
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small and any incorrectness in the theoretical analysis may be a cause for an incorrect 
result. By same arguments one can investigate the effect of the factor T in A/d(V0 on the 
thermodynamics of quasi-2D films. The factor T can be represented as T = T c0 (l + to) 
and one may expect that the usual approximation T « T c0 , which is well justified in 
the Landau theory of phase transitions [2j H] , may be applied. This way of approxima- 
tion can be performed by neglecting terms in the thermodynamic quantities smaller than 
the leading ones. On the other hand practical calculations lead to the conclusion that 
this approximation cannot be made without a preliminary examination because for some 
quasi-2D films it produces a substantial error of about 10% [35]. LT, in which the factor 
T is substituted by T c o is referred to as "simplified Landau expansion" - SLT. All three 
variants of the theory, ET, LT and SLT have been investigated in Ref. [35J. 

2.6 Limitations of the theory 

The general result (1291) - fl32l) for the effective free energy f{ip) has the same domain of 
validity [2] as the GL free energy functional in a zero external magnetic field. When we 
neglect a sub-nano interval of temperatures near the phase transition point we can use 
Eq. (1), provided \to\ — \T — T c0 |/T c0 < 1, or in the particular case of type I supercon- 
ductors, \t Q \ < k 2 [2]. Note, that the latter inequality does not appear in the general GL 
approach. It comes as a condition for the consistency of our approach with the microscopic 
BCS theory for type I superconductors [2]. 

Taking the continuum limit we have to assume that all dimensions of the body, including 
the thickness Lq, are much larger than the characteristic lengths £ and A. The exception of 
this rule is when we consider thin films. Especially for thin films of type I superconductors, 
where ((27r/A) > L 3> a ), we should have in mind that £(T) > A(T), so the inequalities 
£ > A > £o > Ao hold true in the domain of validity of the GL theory |t | < < 1- 
In Ref. [3n] a comprehensive choice of the cutoff A has been made, namely, A = £o (the 
problem for the choice of cutoff A is discussed also in Sec. II. 7 - II. 8). The respective 
conditions for quasi-2D films of type II superconductors are much weaker and are reduced 
to the usual requirements: k > 1/V2, \t \ < 1 and (27r/A) > L ^> a . 

If we do a Landau expansion of i n powers of \ip\ 2 , the condition /) < A 2 should 

be satisfied. In order to evaluate this condition we substitute \ip\ 2 in p = po\ip\ 2 with 
rj 2 = \a\/b, which corresponds to e = 0. As A 2 (T) = 1/p, the condition for the validity of 
the Landau expansion becomes [AA(T)] 2 ^> 1, i. e., (AAo) 2 3> |to|- Choosing the general 
form of A r = (tit/^ ), where r describes the deviation of A T from Ai = A = (7r/£o), we 
obtain (tttk) 2 ^> |£ |; ft — (Ao/£o) is t ne GL parameter. 

Thus we can conclude that in type II superconductors, where k = (Ao/£o) > l/v2, the 
condition (p/A 2 ) <C 1 is satisfied very well for values of the cutoff in the interval between 
A = (7r/£o) and A = (7r/Ao), i.e., for 1 < r < (l/«). For type I superconductors, 
where k < l/y/2 the cutoff value A ~ (l/£o) leads to the BCS condition (|to| < k 2 ) for 
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the validity of the GL approach. Substantially larger cutoffs (A 3> 7r/£ ), for example, 
A ~ (l/^o) f° r type I superconductors with k < 1 lead to a contradiction between the 
BCS condition and the requirement p <C A 2 . 

In our calculations we often use another parameter p T = (1/tctk,) 2 and, in particular, 
fi = fix — (1/ttk) 2 and in terms of p, the condition for the validity of /d^) expansion 
becomes p|to| <C 1, or, more generally, p T |to| ^ 1- Choosing r = 1/ir, we obtain the 
BCS criterion for the validity of the GL free energy for type I superconductors [2]. The 
choice r = (£,o/n\ ) corresponds to the cutoff A T = 1/A . As we will see in Sec. II. 7 
the thermodynamics near the phase transition point in 3D systems has no substantial 
dependence on the value of the cutoff A T but it should be chosen in a way that is consistent 
with the mean-field-like approximation. The cutoff problem for quasi-2D films has been 
investigated in Ref. [35J. It has been shown that the choice A ~ 7r/£ is consistent in the 
framework of GL theory. 

Alternatively, the inequality (p/A 2 ) <C 1 can be investigated with the help of the reduced 
order parameter cp defined by ip = \ip\/r} , where p = p(T = 0) = (aoT^/fr) 1 / 2 is the 
so-called zero-temperature value of order parameter for the GL free energy fo{ip), given 
by Eq. (I3TI) . The reduced order parameter <p will be equal to |£ | for t Q < 0, if only the 
magnetic fluctuations are ignored, i.e., when = rj. Using the notation ip, we obtain 
the condition (p/A 2 ) <C 1 in the form p r <p 2 <C 1. This condition seems to be more precise 
because it takes into account the effect of magnetic fluctuations on the order parameter 



2.7 Bulk superconductors 
2.7.1 Thermodynamics 

The effective free energy f^tp) of bulk (3D-) superconductors is given by Eqs. (1211 - f[3"TT) 
and (148]) . The analytical treatment of this free energy can be done by Landau expansion 
in small (^/~p^\ip\ / h) . Up to order we obtain 

fM « « 3 |^l 2 + |m 4 - qM 3 + f IVf , (so) 



where 



«3 = H TT^T- > ( 51 ) 



h = b+ k ~££, (52) 



2tt 2 A 
6n 



* = ^ , (53) 



and 



_ k B Tpl 

° 3 ~ ~ 6^ • (M) 
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The cutoff A in Eqs. (151 j) - (1541) is not specified and can be written in the form A r = (vrr/^o) 
as suggested in Sec. 2.6. 

It can be shown by both analytical and numerical calculations [31] that |^| 6 -term has no 
substantial effect on the thermodynamics, described by the free energy ( J50l) . That is why 
we ignore this term. The possible phases j^ol are found as a solution of the equation of 
state: 

[dfW/d\1>\ U = • ( 55 ) 

There always exists a normal phase \tpo\ — 0, which is a minimum of fs{ip) for a 3 > 0. 
The possible superconducting phases are given by 




M ± = ^(l±Wl-^j>0. (56) 

Having in mind the existence and stability conditions of iV'oli-phases pf], we obtain that 
the l^oU-phase exists for (16a 3 6 3 ) < 9g 3 and this region of existence always corresponds 
to a minimum of f^{ip). The | ?/>o | --phase exists for < a 3 < 9gf/166 3 and this region of 
existence always corresponds to a maximum of fsiip), i.e., this phase is absolutely unstable. 
For a 3 = 0, |^o|- — an d hence, it coincides with the normal phase. For 9gJ = 16a 3 6 3 
we have |^ |+ = \ip \- = 3g 3 /46 3 and f 3 (\ipo\+) = /aGV'ol-) = 27g|/5126|. Furthermore 
/s(lV'ol-) > for all allowed values of iV'ol- > 0, whereas /s(|'0o|+) < for a 3 < ql/2b 3 , 
and / 3 (|^o|+) > for (g|/26 3 ) < a 3 < 9qf/165 3 . The equilibrium temperature T eq of 
the first order phase transition is defined by the equation f(\il)o\+) — 0, which gives the 
following result: 

26 3 (T eq )a 3 (T eq ) = q 2 3 (T eq ) . (57) 

These results are confirmed by numerical calculations of the effective free energy (1501) [31] . 

The equilibrium entropy jump is AS = VAs and As = — (df 3 (\i/j\) / dT) can be calculated 
with the help of Eq. flBUl) and the equation of state (1551) : 

A S = -|^ | 2 $(|^ |), (58) 



where $(|?/>o|) is the following function: 



\ f , k B Ap . Po k B ,{k B po\ 2 ( rn\ 

The specific heat capacity per unit volume AC = T(dAs/dT) is obtained from (!58j) : 

The quantities As(T) and AC(T) can be evaluated at the equilibrium phase transition 
point T eq , which is found from Eq. (IBTI) : 

T cq 1 fc BPo A , (Po /2 ^/6^) 2 / Tc0 



T c0 ~ 1 2n*a + b + (p^ B /27T 2 A) T c0 \ a ' ' ° ' > 
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provided |AT C | = |T c0 — T eq | -C T c0 . Further we will see that the condition |AT C | <C T c0 
is valid in real substances. The second term in r.h.s. of Eq. flBTl) is a typical negative 
fluctuation contribution, whereas the positive third term in r.h.s. of the same equality is 
typical of first-order transitions 0j. 

To obtain the jumps As and AC at T eq we have to put the solution l^ol+j found from 
Eq. (RSI in Eqs. (15511 - (ISDI). The result will be: 



As = -# U + H^-ltsSL) *kl (62) 



and 

4a / gf c & 



AC^^-^j , (63) 

where b 3c and g 3c are the parameters b 3 and g 3 at T = T eq . As |AT C | = |T c0 — T oq | <C T c0 
we can set T eq « T c0 in r.h.s. of Eqs. (16"2"j) and (I6"3"l) and obtain g 3c = q 3 (T = T cq ) 53(^0) 
and 63c ~ 6 3 (T c0 ). The latent heat Q = — T eq As per unit volume of the first order phase 
transition at T eq can be calculated from Eq. (1621 . 

Now we shall discuss the ratio 

(AT) eq = -9- (64) 

introduced in Ref. [6] as a measure of the temperature size of the HLM effect - a temper- 
ature interval, where the effect can be observed (in [6] the same ratio is denoted by ATi). 
Keeping only the first terms in Eqs. fl62|) - (1631 . we obtain a result for (AT) eq which is four 
times smaller than that for A7\ [6]. To explain the difference let us mention that Eq. (1631) 
gives the jump AC at the equilibrium phase transition point of the first order phase tran- 
sition, described by |^| 3 term [4], while the specific heat jump considered in Ref. [6] is 
equal to the specific heat jump at the standard second order transition AC* = (alT c0 /b) 
[1], and is four times smaller. In fact, neglecting the second term in Eq. (63) and having 
in mind that 6 3 = b, see Eq. (66) below, we have that AC = 4AC. Therefore, we obtain 
that (AT)e q is four times smaller than the respective value A7\ in Ref. [6]. This is valid in 
an approximation in which the second term in Eq. (63) is neglected (this approximation 
is justified in Sec. 2.7.2). The approximation of As and AC by the first (leading) terms 
in Eqs. (62) and (63) does not change essentially these quantities provided the cutoff A 
is comprehensively chosen (see also Sec. 2.7.2). Note that in Refs. [6j [7] as well as in 
all preceding papers, the corrections to the leading terms in As and AC have not been 
taken into account. 



2.7.2 Results for Al 

In order to do the numerical estimates we represent the Landau parameters a® and b with 
the help of the zero-temperature coherence length £0 and the zero-temperature critical 
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magnetic field H c0 . The connection between them is given by formulae of the standard 
GL superconductivity theory [2]: £ 2 = (h 2 / 4ma:oT c0 ) an d H 2 = (4iralT 2 /b). The 
expression for the zero-temperature penetration depth A = (hc/2\/2eH c0 ^ ) is obtained 
from the above relation and Ao = (6/ao^coPo) 1,/2 - We will use the experimental values of 
T c0 , H c0 and £o for Al as given in Table 1. The experimental values for T c0 , H c0 and £o 
vary about 10-15%, depending on the method of measurement and the geometry of the 
samples (bulk material or films) [H] but such deviations do not affect the results of our 
numerical investigations. 



Table 1. Values of T c0 , H c0 , f , «, and \i/j \ for W, Al, In [13]. 



Substance 


T c0 (K) 


H c0 (Oe) 


£o (H 


K 


l^ol x 10" 11 


W 


0.015 


1.15 


37 


0.001 


0.69 


Al 


1.19 


99.00 


1.16 


0.010 


2.55 


In 


3.40 


281.5 


0.44 


0.145 


2.0 



The evaluation of the parameters 03 and 63 for Al gives: 

«s = (« T c0 ) [to + 0.972 x HT 4 (1 + t )r] , (65) 

and 

^ = i + H1Z. (66) 

T 

Setting r = 1 corresponds to the cutoff Ai = (tt/£o) (Sec 2.6). For r = (1/k)ai = 10 2 , 
which corresponds to the much higher cutoff A = (71*/ Ao), we have 63 ~ b, i.e., the p 2 , -term 
in 63, given by Eq. (152]) . can be neglected. However, as we can see from Eq. fl66l) . for r = 1 
the same pg-correction in the parameter 6 3 is of order 0.16 and cannot be automatically 
ignored in all calculations, in contrast to the supposition in Refs. 0, [7]. However, the 
more important fluctuation contribution in 3D superconductors comes from the r— term 
in Eq. (|65|) for the parameter 03. This term is of order 10~ 4 for r ~ 1 and this is consistent 
with the condition \to\ < k 2 ~ 10~ 4 , but for r ~ 10 2 , i.e., for A ~ (tt/Ao) ~ 10 6 /im, the 
same r— term is of order 10 2 , which exceeds the temperature interval (T c0 ± 10 -4 ) for the 
validity of BCS condition for Al (Sec. 2.6). 

These results demonstrate that for our theory to be consistent, we must choose the cutoff 
A r = (vrr/^o), where r is not a large number (r — ► 1 -r- 10). To be more specific, we set 
A = Ai = (7r/£o) as suggested in Ref. [30J. 

The temperature shift t eq = to(T eq ) for bulk Al can be estimated with the help of Eq. flBTl) . 
We obtain that this shift is negative and very small: t eq ~ — 10~ 4 . Note, that the second 
term in the r.h.s. of Eq. fl6T|) is of order 10~ 4 , provided A ~ (l/£o) whereas the third term 
in the r.h.s. of the same equality is of order 10 -5 . Once again the change of the cutoff A 
to values much higher than (tt/£o) will take the system outside the temperature interval, 
where the BCS condition for Al is valid. Let us note, that in Ref. [31] the parameter t 
corresponds to our present notation t - But the numerical calculation of the free energy 
function fs{ip) in Ref. [31] was made for the SLT variant of the theory and the shifted 
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parameter (t Q + 0.972 x 10~ 4 ) was incorrectly identified with t and this has led to the 
wrong conclusion for its positiveness at the equilibrium phase transition point T eq . As a 
matter of fact, the shifted parameter (t + 0.972 x 10~ 4 ) is positive at T eq but t eq = to(T eq ) 
is negative, as firstly noted in Ref. [35] . 

Having in mind these remarks, when we evaluate As and AC for bulk Al we can use 
simplified versions of (1621) and ( 1631) which means to consider only the first terms in the 
r.h.s and to take q 3c ~ q 3 and 6 3c ~ b at T c0 . In this way we obtain 

r ers 

Q = -T c0 As = 0.8 x 10- 2 — -2— , (67) 

LK . cur J 



and 



AC = 2.62 x 10 



erg 
cm 3 



(6* 



The results are consistent with an evaluation of AC for Al as a jump (AC = a^Tco/b) at 
the second order superconducting transition point [6] that, as we have mentioned above, 
is four times smaller than the jump AC given by Eq. (1681) . 

A complete numerical evaluation of the function f^{ip) and the jump of the order param- 
eter at T eq for bulk Al was presented for the first time in Ref. [3 lj . The results there 
confirm that the order parameter jump and Q for bulk type I superconductors are very 
small and can hardly be observed in experiments. 

We will finish the presentation of bulk Al with a discussion of the ratio It can be 

also written in the form 

32! fk%1%\ ( e 2 



and it differs by a factor 1/4 from the respective result in Ref. (Sj. The reason for this 
difference was explained at the end of Sec. 2.7.1. From Eq. ([6"9"j) we have 

(AT) C q = 6.7 x 10- 12 (TX 2 £o 6 ), (70) 

and multiplying the number coefficient in the above expression by four we obtain Eq. (10) 
presented in Ref. [6]. 



2.8 Quasi-2D films 

Following Refs. (32 33J we can present the free energy density f(ip) — (F(ip)/LiL 2 ) in 
the form 

f(<p) = ^ [2t ^ 2 + ^ + C(l + t )r(/V)] , (71) 

where 

r(y) = (l+y)ln(l + y)-ylny, (72) 
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To obtain Eqs. (171]) - (172]) we have set A = (7r/£ ) and introduced the notation ip = 
IV'I/W ?7o is defined in Sec. 2.5. Some of the properties of free energy (171]) were analyzed 
in Ref. [32] for Al films and in Ref. [33] for films of Tungsten (W), Indium (In), and 
Aluminium (Al). Here we will briefly discuss the main results. 

The equilibrium order parameter (fo > corresponding to the Meissner phase, can be 
easily obtained from the equation df(cp)/dcp = and Eq. ( ITT]) : 

^WUWJ, + M . (73) 

2 V Wo/ 

The logarithmic divergence in Eq. (173]) has no chance to occur because <y?o is always 
positive and does not tend to zero. 

The largest terms in the entropy jump 5 s and the specific heat jump 5C at the equilibrium 
first order phase transition point T eq are given by [321 E3] 

'—si/- (74) 

where ^ eq = <fo(T = T eq ), and 

= . (75) 

4ttT c0 

The latent heat of the phase transition [1] is given by Q = — T cq 5s and Eqs. (]73"]) - (j7I]) . 
Since the temperatures T eq and T c0 have very close values, the difference between the 
values of Q, 5s, and 5C at T c0 and T eq , respectively, can also be ignored, for example, 
\5C(T cq ) - 5C(T c0 )\/5C(T c0 ) < 1 and we can use either 5C(T c0 ) or 5C(T cq ) [32]. 

The equations (171]) - (1731) corresponding to quasi-2D films are quite different from the 
respective equations for bulk (3D-) superconductors but it is easily seen that the relatively 
large value of the order parameter jump (p 2 in thin films again corresponds to relatively 
small value of the GL parameter k. That is why we consider element superconductors 
with small values of k and study the effect of this parameter, the critical magnetic field 
H c0 and the film thickness L on the properties of the fluctuation-induced first order phase 
transition. 



We use experimental data for T c0 , H cQ , £ and k for W, Al, and In, published in Ref. 
(see Table 1). In some cases the GL parameter k can be calculated with the help of 
the relation k = (Aq/£o) and the available data for £ and Ao- In other cases it is more 
convenient to use the following representation of the zero-temperature penetration depth: 

Ao = • (76) 

2 v / 2e# c0 £o 

The value of |"0o| i n Table 1 is found from 

/ rn \ 1/2 

m = hw ^ co - (77) 
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Figure 2: Order parameter profile <p(to) of Al films of different thicknesses: L = 0.05 /J,m 
("+"-line), L = 0.1 fim (o), and L = 0.3 //m (•) [33]. 

The order parameter dependence on the reduced temperature difference to is shown in 
Fig. 2 for Al films of different thicknesses. It is readily seen that the behavior of the 
function <^o(^o) corresponds to a well established phase transition of first order. The 
vertical dashed lines in Fig. 2 indicate the respective values of t eq = to(T cq ), at which 
the equilibrium phase transition occurs, as well as the equilibrium jump <p (T eq ) = <p eq 
for different thicknesses of the film. The parts of the y9o(to)-curves, which extend up to 
to > t eq , describe the metastable (overheated) Meissner states, which can appear under 
certain experimental circumstances (see in Fig. 2 the parts of the curves on the r.h.s. of 
the dashed lines). The value of (p eq and the metastable region decrease with the increase 
of the film thickness, which shows that the first order of the phase transition is better 
pronounced in thinner films and this confirms a conclusion in Ref. [32] . 

These results are justified by the behavior of the free energy as a function of to- We used 
Eqs. (!7T|) - (173|) for the calculation of the equilibrium free energy f[fo(to)]. The free energy 
for Al films with different thicknesses is shown in Fig. 3. The equilibrium points T eq of the 
phase transition correspond to the intersection of the /(</?o)-curves with the to-axis. It is 
obvious from Fig. 3 that the temperature domain of overheated Meissner states decreases 
with the increase of the thickness L . 

The shape of the equilibrium order parameter y?o(^o) i n a broad vicinity of the equilibrium 
phase transition for thin films (Lq = 0.05/xm) of W, Al, and In was found from Eq. (1731) . 
The result is shown in Fig. 4. The vertical dashed lines in Fig. 4 again indicate the 
respective values of t eq = t (T eq ), at which the equilibrium phase transition occurs as well 
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Figure 3: The free energy /(to) for Al films of thickness: Lq = 0.05 fxm ("+"-line), 
L = 0.1 /im (o), L = 0.3 /im (•) [33]. 
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Figure 4: Order parameter profile <£>(to) of films of thickness Lq = 0.05 /im: W ("+"-line), 
Al (o), and In (•) [33]. 
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as the equilibrium jump ip (T eq ) = <p eq in the different superconductors. 

The order parameter jump at the phase transition point of In (the In curve is marked by 
points in Fig. 4) is relatively smaller than for W, and Al, where the GL parameter has 
much lower values. The same is valid for the metastability domains; see the parts of the 
curves in Fig. 4 on the left of the vertical dashed lines. It is obvious from Fig. 4 and Table 2 
that the equilibrium jump of the reduced order parameter (p eq of W has a slightly smaller 
value than that of Al, although the GL number k for W has a ten times lower value, 
compared to k of Al. Note, that in Fig. 4 we show the jump of <£> eq , but the important 
quantity is \?p\eq = iV'olv^eq- Using the data for L = 0.05/xm from Tables 1 and 2 we find 
for |^| eq the following values: 0.1 x 10 11 for Al, 0.05 x 10 11 for In, and 0.02 x 10 11 for W. 
This result shows that the value of the critical filed H c q is also important and should be 
taken into account together with the smallness of GL number when the maximal values 
of the order parameter jump are looked for. Thus the value of the order parameter jump 
at the fluctuation- induced phase transition is maximal, provided small values of the GL 
parameter k are combined with relatively large values of the critical field H c q. In our case 
Al has the optimal values of these two parameters. 



Table 2. Values of t eq , f e q, and Q (erg/cm 3 ) for films of W, Al, and In with different 
thicknesses Lq (/jm) [33J. 







Al 






In 






w 








teq 




Q 


teq 


Veq 


Q 


teq 


</>eq 


Q 




0.05 


-0.00230 


0.041 


1.95 


-0.00167 


0.025 


3.94 


-0.00174 


0.039 


1.6 x 10" 


-4 


0.1 


-0.00147 


0.032 


0.80 


-0.00094 


0.017 


1.82 


-0.00118 


0.032 


1.1 x 10" 


-4 


0.3 


-0.00070 


0.023 


0.41 


-0.00037 


0.010 


0.63 


-0.00064 


0.023 


5.6 x 10" 


-5 


0.5 


-0.00048 


0.016 


0.20 


-0.00029 


0.008 


0.40 


-0.00048 


0.020 


4.1 x 10" 


-5 


1 


-0.00029 


0.012 


0.11 


-0.00013 


0.006 


0.23 


-0.00032 


0.016 


2.7 x 10" 


-5 


2 


-0.00017 


0.009 


0.06 


-0.00008 


0.004 


0.10 


-0.00021 


0.013 


1.8 x 10" 


-5 



The importance of the zero-temperature critical magnetic field H c0 for the enhancement 
of the jumps of the certain thermodynamic quantities at the equilibrium phase transition 
point T eq becomes obvious from Eqs. (1741 . (1751) and (1771) . Eq. (ITTj) shows that the order 
parameter jump |-0| eq = |^ |<^ eq is large for large values of H c0 and £ - This is consistent 
with the requirement for relatively small values of the GL parameter k, as shown by 
Eq. (1761) . Therefore, the unmeasurable ratio Q/5C discussed in Ref. [6] does not depend 
on the value of the critical field H c o but the quantities Q and 5C themselves, as well as the 
order parameter jump \ip\ eq depend essentially on H c0 . The values of the reduced order 
parameter jump <p eq for films of Al, In and W of the same thickness have the same order 
of magnitude while the respective order parameter jump \ip\ eq is one order of magnitude 
higher for Al than for W, as shown above. The effect of the critical magnetic field H c q on 
the latent heat Q is, however, much stronger and, as is seen from Table 2, the latent heat 
Q in W films is very small and can be neglected while in Al and In films it reaches values, 
which could be measured in suitable experiments. This is so because the latent heat is 
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proportional to the difference [H 2 Q /8n ~ ^l^ol 4 ] between the energies of the ground state 
(superconducting phase at T = 0) and the normal state. It is consistent with the fact that 
the fluctuation contribution to the free energy, i.e., the C— term in the r.h.s. of Eq. flTTj) is 
generated by the term of type \ip\ 2 J d D xA 2 (x) in the GL free energy. At T = this free 
energy term is also proportional to the mentioned difference between the free energies of 
the superconducting ground state and the normal state. 

The shift of the phase transition temperature t eq = |(T cq — T c0 )|/T c0 , the reduced value (p eq 
of the equilibrium order parameter jump \if)\ eq , and the latent heat Q of the equilibrium 
transition are given for films of different thicknesses and substances in Table 2. This table 
shows that the shift of the phase transition temperature is very small and can be neglected 
in all calculations and experiments based on them. The values for ip eq for different Lq 
and those for |f/> | given in Table 1 confirm the conclusion, which we have made for films 
of Al, In, and W with L = 0.05/im. The latent heat Q has maximal values for In, where 
the critical field is the highest for the considered materials. 

2.9 Final remarks 

Our analysis shows that the MF studies of the HLM effect have a well defined domain of 
validity for both 3D- and quasi-2D superconductors. Our conclusion is that the MF theory 
of the magnetic fluctuations in superconductors and, in particular, the MF prediction for 
the fluctuation driven weakly first order phase transition in a zero external magnetic field 
in bulk superconductors [6] and quasi-2D superconducting films j3(2, [32], [33] is reliable and 
can be tested by experiments. While the HLM effect in bulk systems is inobservantly 
small, in quasi-two dimensional superconductors this effect is much stronger and may be 
observed with available experimental techniques. 

Our consideration of quasi-2D superconductors is highly nontrivial in view of the rele- 
vance of the effective Landau parameters dependence on the thickness of the films, Lq. 
We have justified this dependence by simple heuristic arguments and by a reliable con- 
sideration of the 2D-3D crossover. In contrast to initial expectations [6] that films made 
of superconductors with extremely small GL parameter k such as Al and, in particular, 
W will be the best candidates for an experimental search of the HLM effect, our careful 
analysis (firstly published in Ref. [33J ) definitely gives somewhat different answer. The 
Al films still remain a good candidate for transport experiments, through which the jump 
of the order parameter at the phase transition point could be measured, but surprisingly 
the W films turn out inconvenient for the same purpose, due to their very low critical 
field H c q. The importance of the critical magnetic field H c q for the clearly manifested 
first-order phase transition has been established and discussed, too. Although In has ten 
times higher GL number k than Al, the In films can be used on equal footing with the 
Al films in experiments intended to prove the order parameter jump. Here the choice of 
one of these materials may depend on other features of experimental convenience. As far 
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as caloric experiments are concerned, the In films seem to be the best candidate due to 
their high latent heat. 

As shown in this review experiments, designed to search the HLM effect, can be performed 
by both type I and type II superconductors. In experiments the external magnetic field 
cannot be completely eliminated. Then vortex states may occur for H = \H\ > below 
T c = T C (H) < T c0 in type II superconducting films and this will obscure the HLM effect. 
Note, that in both type I and type II superconductors the external magnetic field H 
generates additional entropy jump at the phase transition point T C (H) and this effect 
can hardly be separated from the entropy jump ( |74l) caused by the magnetic fluctuations 
in the close vicinity of T c q. Therefore, in experiments intended for searching of HLM 
effect, the external magnetic field should be minimized as much as possible. For quasi- 
2D superconductors, where the HLM effect is relatively strong and the latent heat can 
exceed several ergs, one should ensure external fields less than 10 Oe, or, in more reliable 
experiments, less than 1 Oe |35j . 

The temperature range around T c q, where the HLM effect is significant, has been estimated 
in Ref. [6] to a few microdegrees, and later, this estimate has been confirmed in a more 
accurate way [15]. According to our point of view, the possibilities for the observation 
of the effect are related mainly with its magnitude, because the range of temperatures 
where it occurs always exists, even in the critical region of strong ^-fluctuations, as the 
RG studies indicate (Sec. III). The review of the results for Al, In, and W in Sec. II shows 
that the metastability domains (of overheating and overcooling) are much larger than the 
Ginzburg critical region. This result justifies the reliability of the MF treatment. As 
one may see from Figs. 2-4 the metastability temperature interval is relatively larger for 
smaller values of the GL number k, but in order to ensure large values of some measurable 
thermodynamic quantities as, for example, the latent heat Q and the specific heat capacity 
C we must choose a material with large critical field H c . The experimental verification 
of the order parameter jump can be made by transport experiments. As we see from 
Figs. 2-4 and Table 2, this quantity has maximal values for W, where the parameter k is 
very small. Having in mind also the large metastability regions in this material, one may 
conclude that W is a good candidate for a testing the HLM effect by transport experiments 
(measurements of the superconducting currents), provided some specific disadvantages of 
this material (quite low T c q, etc.) from the experimental point of view do not contradict 
to this our conclusion. Another effect, which is relevant to the present discussion, is the 
known variation of the GL parameter k with the variation of the thickness L (see, e.g., 
Ref. [16] . The parameter k < l/y/2 of a type-I bulk superconductor may change up to 
values corresponding to a type-II superconductor with the decrease of the thickness L 
below 10 _7 m. 

In Ref. [47] a very interesting interrelationship between the first order phase transition 
and the vortex fluctuations in type-I 2D superconductors has been investigated by both 
analytical and numerical methods. A recent study [H] of HLM effect in superconducting 
films has been proven wrong in our comment |49j; see, also, the Corrigendum [50] to the 
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work [48 J . In further papers the same authors and co-workers confirm [5TI 152"] our results, 
published in Ref. [301 Ell E2l [33J, [MJ [35] , and discussed in this Section. Another paper 
has been intended to the treatment of gauge effects in superconductors with the help of a 
Gaussian effective functional [53], which gives results quite near to mean-field ones, and 
is wrong within the present discussion. 

3 RG STUDIES AND OTHER TOPICS 

3.1 Notes about RG 

The RG methods are intended to investigation the strong fluctuation interactions in the 
critical region of second order phase transitions and multi-critical phenomena [HE]. In 
studies of complex systems with several orderings or with influence of additional gauge 
fields, as in Eq. (1), usually the k— space RG methods are used because of their wider 
applicability. The methods differ from each other in specific technical details but all of 
them lead to the same predictions about the phase transition properties. The common 
feature is that the mentioned methods are based on the so-called loop expansion [H, [5J, [12] . 
The RG equations are, in fact, infinite asymptotic series, which are truncated at one-, 
two-, or, in rare cases, higher orders in the loop expansion. The total infinite series 
can be summed only in case of trivial (usually exactly solvable) models and in some very 
exceptional cases as, for example, the RG series for interacting real bosons in the quantum 
limit T -> [53]. 

For the simple 4 — theory the RG series can be derived and analyzed to high orders in 
the loop expansion whereas for more complex models this can be practically done within 
the one- and two-loop approximations. As the RG series are asymptotic, normally, the 
one- and two- loop orders give all important features of the specific system of interest: (i) 
the presence of stable fixed points (FP) of RG equations and the associated with them 
types of critical behavior, and (ii) the lack of stable FPs and related conclusions about 
the lack of standard critical or multi-critical behavior. The higher orders of the theory 
are more relevant to investigations of the asymptotic properties of the loop series than for 
obtaining of qualitatively new critical behavior, or, other new qualitative characteristics of 
the system. The latter are reliably obtained within the one- and two-loop approximations. 

In our consideration we shall give examples of RG equations in the original variant of 
the Wilson-Fisher RG approach; see, e.g., Ref. [4]. It is based on the loop expansion 
and an explicit scaling procedure for the wave vectors (k = e~ l k') and the fields ip(k) = 
e (i-W 2 ) z ^/(£') an d Aj(k) = e { - l ~ VA ^ 2 ' )l A'^k'), where e l > 1 is the so-called re-scaling factor 
(/ > 0), and rjip and t\a are the anomalous dimensions (alias Fisher exponents) of the fields 
ip(k) and Aj(k), respectively. 
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3.2 The order of the phase transition 



In Ref. [6] the simultaneous effect of if)- and A-fluctuations in the fluctuation Hamiltonian 
(1) has been investigated by one-loop RG and e = (4 — <i)-expansion [4] for the general 
case of n/2-component complex vector field ip = {ip a ]& = 1, w/2}; note, that this field 
is equivalent to a n-component real vector field. It has been shown [6] that a stable FP 
exists below four dimensions d — (4 — e) only for symmetry indices n > 365.9, which are 
far above the real symmetry index n = 2 for usual superconductors and numbers n = 4, 6 
corresponding to superconductors with certain unconventional Cooper pairings [55J. For 
n < 365.9 real FP does not exist at all. Besides, the 1/n-expansion has been used for 
the calculation [6] of the critical exponents in the so-called large- n limit (alias "Hartree 
limit") [HE], where (n > 365.9) the phase transition is definitely of second order. 

According to the normal interpretation of RG results, the lack of stable FP is an indication 
for a lack of standard second order phase transition. But usually, additional (non-RG) 
arguments are needed to determine the actual order of the phase transition. In the 
present case, the lack of stable FP could be a result of the same mechanism that produces 
a fluctuation- driven weakly first order phase transition in MF approximation (Sec. 2). As 
RG takes into account both fluctuations of ip and A, one can conclude that the result 
for the weakly-first order phase transition in a zero external magnetic field is valid for 
both type I and II superconductors. Of course, the size of the effect (the size of jumps of 
energy, latent heat, order parameter) will depend on the specific substance. This path of 
investigations has been further developed in the paper [7J and several problems, opened 
in the short Ref. [6], have been solved. 

In Ref. [40], where both 2D and 3D superconductors are considered, the HLM effect 
has been confirmed in one-loop order (annealed disorder [40J) as well as the availability 
of second order phase transition for large-n has been proven in another variant of the 
Hartree limit (n — > oo). 

Next, the effect has been confirmed in a RG investigation of the equation of state below 
the phase transition point T c0 [56J, [57]. These investigations show that the effect may 
occur under quite restricted conditions for the vertex parameters e and b of the functional 
(1), in particular, when the effective charge \e*\ exceeds some value. A strong restriction 
on the HLM effect, but in the opposite direction - a requirement for a sufficiently small 
effective electric charge |e*|, has been found also by a Monte Carlo (MC) study [58J of a 
lattice version of the model (1). Another MC study [59] concludes that at fixed effective 
charge (|e 2 | = 5) the HLM effect strongly depends on the GL parameter k = A/£: it is well 
established for k <C 1, becomes very weak for n ~ 1 and vanishes for large k. On the basis 
of this picture, a proposal is made [59J that a tricritical point exists at some "tricritical" 
value of k. As k is related with the charge e, this proposal seems to be in conformity with 
other investigations. Note, that the first prediction for change of the phase transition 
order with the variation of k and for the possibility of "tricritical point " to occur at some 
value of k, has been made in Ref. [SU] on the basis of duality arguments |2] and analytical 
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calculations; see further development of this approach in Refs. [621 E21 EH E2] • A further 
MC study [66] of type II superconductors indicates that these systems undergo a second 
order phase transition of universality class 3D XY rather than a weakly first order phase 
transition. The two-loop RG studies [671 EH] show that the account of higher orders of the 
loop expansion allows a lowering of the critical value n c = 365.9, but the RG equations 
still have no real FP at D = 3 and n = 2 (conventional 3D superconductors), except for 
a treatment by Pade analysis [68] (see, also, the review |37J). In 2 + e dimensions, the 
nonlinear a model exhibits a second-order phase transition for all values of n > [69] . 
Owing to this result and more extensive studies [ZD, [72] of the dependence of the effective 
free energy on the GL parameter k, one may suppose that the critical value n c vanishes 
at some dimension 2 < d < 4, that is, that a tricritical point "located" at some specific 
value of k exist where the superconducting phase transition changes its order; see also, 
Ref . [70] . RG studies [721 E3] a t fixed dimension d have also been carried out with the 
aim to determine the phase transition order (see, also. Ref. [70] for a comment). 

3.3 Disorder effects 
3.3.1 Quenched impurities 

Here we will discuss mainly disorder described by quenched impurities and inhomo- 
geneities with Gaussian distributions (see, e.g., Ref. [1]). The effects of annealed and 
quenched disorder in classical versions of Abelian-Higgs models, equivalent to the GL 
model (1) have been investigated by Hertz [7H [75] in the context of the theory of spin 
glasses in case of Dzyaloshinskii-Moriya interaction. The specific feature of this approach 
is that the disorder is associated with the vector gauge field A(x) and can be used to 
describe superconductors in random, uncorrelated, or, which is the same, short-range (S-) 
correlated magnetic fields rather than usual ones. In high energy physics this approach 
makes contact with gauge fields, where a Higgs field is coupled to a random color field. It 
is natural to expect, as has been proven [101 El], that the annealed gauge model will bring 
a fluctuation-driven first order phase transition at usual symmetry indices and a second 
order phase transition in the Hartree limit (n — ► oo). In Ref. [10] the phase transition 
in case of quenched impurities has been predicted to be of second order in a calculation 
to one-loop order, in contradiction to the conclusion for a lack of stable nontrivial FP of 
the RG equations within the same one- loop order [71]. The origin of this discrepancy is 
clear from the argument that the lack of stable FP may result without any change of the 
Hamiltonian structure, whereas the conclusion for the second order of phase transition 
in Ref. [10] has been made, perhaps, incorrectly, based on the argument of the |?/>| 4 -term 
absence for D = 3 and the |?/>| 2 ln|?/>|-terni absence for D = 2. 

More usual case of quenched disorder in Abelian-Higgs models has been considered in 
Refs. [76"! [77] . This is the case when the parameter a in (1) has a random part <p(x), 
i.e. a is substituted with a = a + ip(x), where the Fourier amplitudes S(k) of the random 
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function <p(x) obey, for example, the Gaussian distribution 



(<p(k)<p(k% and = A(k)S- 



fc,fe' 



(78 



Here the parameter A(k) > is a non-negative function for any < k < A (hereafter, 
A = 1). This disorder is of type random impurities, or, what is the same, "random critical 
temperature" [I]. In case of the so-called "uncorrected", or, (5-) short range correlated 
quenched impurities [A(k) = A = const], the system exhibits a spectacular competition 
between the impurities and gauge effects. We shall show this by the example of the RG 
equations derived in Ref. [77]. Using units, in which kg = T c0 = A = 1 and the notations 
go = (2e/hc) and 7 = h 2 /4m the RG equations corresponding to the impurity model given 
by Eqs. (1) and (78) can be written in the form 
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When the disorder is neglected (A = 0; pure superconductors) the Eqs. (79) - (83) 
coincide with those for pure superconductors |6j and give the results briefly discussed in 
Sec. 3.2. For impure superconductors (A > 0), these equations describe a new stable 
FP of focal type [77] for dimensions D < 4, which exists for symmetry indices n > 1 
and has a physical meaning at dimensions D > D c (n) = 2{n + 36)/(n + 18) [77J. In 
the impure superconductor, the new focal FP occurs exactly in the domain of symmetry 
indices n < 365.9, where the HLM analysis [6] yields a lack of FP for the respective 
pure system (for details, see Ref. [77]). Having in mind the asymptotic nature of the e- 
expansions within the RG approach [H [78] , one may conclude that this focal FP governs 
the critical behavior at the real dimension D = 3 in conventional superconductors (n = 2), 
although the direct substitution of n = 2 in the above expression for the "lower critical" 
dimension D c (n) yields D c (2) = 3.8. This problem has been further investigated in 
Ref. [79J. Long-range quenched disorder within the same RG approach to the model (1) 
has been considered in Ref. 



3.3.2 Random field disorder 

The effect of another relevant type of quenched disorder - the so-called "random field 
disorder" [4] - on the phase transitions described by the GL functional (1) has been 
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considered in Ref. [SI]. In order to describe this type of disorder one should add to the 
GL model (1) the following terms |81j : 

_/^W)+c.c.], (84) 

where h = {h a ; a = 1, . . . , n/2} is a random field with a Gaussian distribution of the type 

(h* a (k)hp(k')) rand = d^ghlk- 9 . (85) 

In Eq. (85), 8 = describes the k ~ asymptote of short-range random correlations (in x- 
space, \x—x'\ a ~ d ; a < 0), and < 9 < d describes long-range correlations (0 < 6 = a < d); 
see Refs. [H [82], [83] . 

In both short-range and long-range random field correlations, a new stable FP has been 
obtained [8~T] . For short-range random correlations this FP is stable below the upper 
critical dimension Dy = 6 and for symmetry indices n > n c = 10. This FP describes a 
quite specific critical behavior. The situation for long range random correlations is more 
complicated and here we will advise the reader to follow the discussion in Ref. [81J. 

3.3.3 General conclusion 

The annealed disorder leaves the phase transition properties almost identical to those in 
the pure Abelian-Higgs model (1), whereas the quenched disorder gives a lower critical 
values of the symmetry index n, below which the weakly first order phase transition (HLM 
effect) exists. Thus one may conclude that the quenched disorder acts in a direction of 
"smearing" the fluctuation-driven first order phase transition, i.e. the HLM effect seems 
to be much weaker in such systems. 

3.4 Liquid crystals 

The first theoretical study of the HLM effect on the properties of the nematic-smectic 
A phase transition in liquid crystals has been performed in Refs. [22l [23] . The main 
result of these investigations is the confirmation of the HLM effect in a close vicinity of 
this liquid crystal phase transition. Perhaps, because of the liquid crystal anisotropy, 
which explicitly enters in the propagator of the gauge field [23J, here the critical value 
n c = 38.17 of the symmetry index n, below which the weakly first order transition occurs, 
is lower than in pure superconductors {n c = 364.9) but it is still quite above the real value 
n = 2. Thus the RG predictions in Refs. [221 E3] in seventies of the previous century 
were in favor of a weakly first order phase transition whereas the experimentalists [SH 
[85] at the same time claimed that the same nematic-smectic A phase transition is of 
second order (for triple and tricritical points in nematic-smectic- A-smectic-C systems, and 
the interrelationship to the present problem, see Ref. [86]), as well as later papers [871 
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[55] . While in 3D superconductors the HLM effect is very weak and unobservable with 
available experimental techniques, the size of the same effect in 3D liquid crystals allows an 
observation in case of very precise experiment and elimination of obscuring effects as, for 
example, neighborhood to tricriticality and liquid crystal anisotropy. Experiments [89], 
190] consistent with the HLM effect has been carried out but, as noticed in Ref [91], 
these experiments are at their resolution limits and the full implications of the HLM 
theory remains to be tested. On the other hand, experiments in Ref. [91] demonstrate a 
discrepancy with the simple approximation (ip =const) for the smectic-A filed within the 
mean-field like approach (Sec. II). 



3.5 Quantum effects 

The first account of quantum correlations ( "fluctuations" ) [H [92} 93J on the properties of 
the phase transition to superconducting state has been performed in Ref. [94J . It has been 
shown [9l] that the dynamical critical exponent z, produced by the quantum correlations 
at finite temperatures (T > 0), is given by z = 2 + 18e/n. This is purely gauge result 
because, as shown in Ref. [95], the critical dynamics of a superconductor when neglecting 
the local gauge effects (A = 0) is identical to that described by the time dependent GL 
model (TDGL) in case of lack of conservation laws [96]; a prediction of the latter result 
has been given for the first time in Ref. [92] . 

Another paper [97] on a gauge model of type (1) is intended to the investigation of the 
quantum phase transition (T c0 = 0) in granular superconductors. The simultaneous effect 
of the local U(l) symmetry, disorder and quantum fluctuations has been considered in 
the Ref. [98, 99J, where disordered electronic systems at T = are studied. The problem 
was further extended to the treatment of quantum phase transitions in underdoped high- 
temperature superconductors [100[ 1101] . 

In Ref. |102j the quantum phase transition at zero temperature from superconducting 
state to the state of normal metal has been discussed within the framework of the HLM 
approximation (Sec. 2.2). In the limit T — > 0, a dimensional crossover of type d — > (d + 1) 
occurs due to the quantum fluctuations of the gauge field A(r, x), where r is the imaginary 
time [U [93]. Thus the effective free energy of the three dimensional "quantum system" 
at T = corresponds to the free energy of the respective four- dimensional (4D) classical 
system. Remember, that the latter is described by Eqs. (29)-(31) and (49). Thus the 
MF like analysis of 3D superconductor in the zero-temperature limit T — > will give a 
fluctuation-driven first order phase transition, described by a logarithmic term of type 
j^^lnl^l rather than the |-?/>| 3 -term known from Ref. [6]. Accordingly, in the same zero- 
temperature limit (T c = 0) the 2D- and quasi-2D superconductors will be described by an 
effective free energy, containing a |^| 3 -term as is in 3D superconductors with T c ^ 0. The 
complete analysis of the quantum phase transitions (T c = 0) of type-II superconductors 
beyond the MF-like approximation (Sec. 2.2) requires the account of both thermal and 
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quantum fluctuations of the order parameter field ip{r, x) with the help of RG; for example, 
such fluctuations have been systematically taken into account in Refs. [9?] 198]. 

3.6 Complex systems 

The RG investigation [103J of systems with superconducting and other (non-magnetic) 
orderings shows that the gauge field A(x) leads to a drastic modification of the critical 
behavior in a close vicinity of bicritical and tetracritical points (for such points, see, e.g., 
Ref. OH])- Near these multicritical points the superconducting fluctuations can be en- 
hanced by the fluctuations of another ordering field and, hence, the HLM effect in such 
complex systems seems to be stronger than in the more simple case, discussed so far 
|104j . The coupling between the superconducting field ip and the magnetization mode 
M(x) in models of ferromagnetic superconductors intended to describe the coexistence 
of superconductivity and ferromagnetism in ternary compounds, has been investigated 
by RG in Ref. [1051 1106] . It has been shown again that the gauge field A(x) produces a 
weakly first order phase transition, enhanced by the M-fluctuations. The same problem 
has been further opened in Ref. [107] by a lattice version of the GL model and dual argu- 
ments [61] with the conclusion that under certain circumstances the phase transition from 
the disordered phase to the phase of coexistence of ferromagnetism and superconductivity 
should be, at least in type II superconductors, of second order in contrast to the RG pre- 
diction [1051 1106] . Another investigation of coupled magnetic and superconducting order 
parameters, where the gauge effects are also relevant, is intended to describe a state of 
vortex solid in rear-earth and borocarbide compounds [1081 1109] . 

3.7 Unconventional superconductors 

Another interesting problem is the behavior of unconventional superconductors [551 1110] , 
where the field ip is a complex vector: ip = {ip a , a = 2, 3}, i.e., n = 4 or 6, depending on 
the type of the unconventional Cooper pairing. As already mentioned, in RG studies we 
may consider an arbitrary value of the symmetry index n/2 > 1, and to put n/2 = 2,3 
in our discussion of real unconventional superconductors. In order to treat these systems 
we must add two new terms to the free energy (1), namely, 

n/2 

^(^T + ^£hm (86) 

where the parameter u describes the anisotropy of the unconventional (non-s-wave) Cooper 
pairing and the w-term has been included to represent the crystal anisotropy; as given in 
(86) this term represents a cubic crystal anisotropy. It has been shown [lllj, that HLM 
FP in this case is unstable towards the parameters u and v, describing the anisotropy 
of the Cooper pair and the crystal anisotropy, and the new FP points that appear in 
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the theory are unstable for usual values of the symmetry index (n = 4,6). Analytical 
calculations [111] have been used to demonstrate that one of the FP points that appear 
in the RG equations, exists for n > n c = 10988 and is stable for some values above n c , 
at least, in the Hartree limit (n — > oo). Owing to these results, and mainly due to the 
very high value of n c , a conclusion has been made that the HLM effect should be more 
pronounced in unconventional superconductors [lllj . 

A next step along this path of RG studies has been performed in Ref. |112j . where un- 
conventional superconductors with quenched impurities have been investigated by RG to 
one-loop order. For the most interesting cases of n = 4 and n = 6, the RG equations 
have a focal FP, which is known from the results for conventional (u = v = 0) impure 
superconductors (Sec. 3.3.1). However, here this FP is unstable towards the parameters 
u and v, which means that the superconducting phase transition in a close (asymptotic) 
vicinity of the phase transition point is not a conventional second order phase transition 
but rather this transition is a modified version of the fluctuation-driven first order phase 
transition, known from Ref. [6]. 

Another interesting problem is related with the gauge effects in ferromagnetic unconven- 
tional superconductors such as UGe2, URhGe, ZrZn 2 , where a rich picture of various phase 
transitions and (multi) critical points occurs [113] . including the existence of a completely 
new type of finite temperature critical behavior [1141 I115[ I116j and quantum critical 
points. This problem is quite new and since now there are no investigations of the effects 
of the gauge field fluctuations 5A(x) on the phase transition properties in ferromagnetic 
unconventional superconductors. 

3.8 External field effect 

In real experiments the external magnetic field can hardly be completely eliminated. The 
study (38] of the HLM effect in the wider context of nonzero external magnetic field 
H , shows that the weakly first phase transition, discussed so far can be obtained in 
Landau gauge from a renormalized theory of the equation of state in one loop order; see 
also Ref. [45]. The phase transition at the second critical magnetic field has been 
studied by a e = (6 — d) expansion within the one loop approximation with the conclusion 
that the magnetic fluctuation effects should produce a fluctuation induced first order 
phase transition [117j . Further investigations of this problem have been performed in the 
Hartree limit {n — > 0) [118[ I119j . and in higher orders of the loop expansion |12U[ 1121] ; 
for comments, see Refs. |122[ 1123] . 
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